Stripes, topological order, and deconfinement in a planar t-Jz model 
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We determine the quantum piiase diagram of a two-dimensional bosonic t-Jz model as a function of the lattice 
anisotropy 7, using a quantum Monte Carlo loop algorithm. We show analytically that the low-energy sectors 
of the bosonic and the fermionic t-Jz models become equivalent in the limit of small 7. In this limit, the ground 
state represents a static stripe phase characterized by a non-zero value of a topological order parameter. This 
phase remains up to intermediate values of 7, where there is a quantum phase transition to a phase-segregated 
state or a homogeneous superfluid with dynamic stripe fluctuations depending on the ratio Jz/t. 

PACS numbers: 05.30.-d, 05.70.Fh, 05.30.Jp, 75.10.Jm 



0) 



I 

o 
o 



> 

O 
m 
o 



o 
o 



X 



During last decade a lot of attention has focused on the 
study of inhomogeneous structures in strongly correlated ma- 
terials, such as the copper oxide based high-temperature su- 
perconductors (HTSC) \M- All HTSCs share a common fea- 
ture - an antiferromagnetic (AF) insulating parent state which 
evolves into a variety of phases upon doping with carriers, 
either chemically or using some external probe. In partic- 
ular, there is experimental evidence fS] that in part of the 
phase diagram some of these compounds feature inhomoge- 
neous charge and spin textures, commonly known as stripes. 
Existence of such textures is usually justified by competing 
interactions between the particle constituents, which lead to 
new locally phase-separated states, all of them characterized 
by order parameters (OPs) implying that a broken (e.g., trans- 
lational) symmetry state is in place. The real interest lies in 
the relation between these stripe phases and superconductiv- 
ity, since it is unknown how these quantum orders cooperate 
with or compete against each other |31 . 

In the present manuscript we argue for the existence of 
static and dynamic stripe phases in certain (2H-l)-dimensional 
lattice models of strongly correlated materials with no charge- 
ordered state associated to it, contrary to cuiTent understand- 
ing. Our starting point is not an assumed system of interacting 
stripes but a system where the latter emerge from the com- 
petition between antiferromagnetism and delocalization. We 
demonstrate that the essential physics is related to the exis- 
tence of spin antiphase domain structures which are known 
to be ubiquitous in various doped AF insulators. The t-J^ 
Hamiltonian is a minimal model for studying that competition 
which is also present in the t-J model. Particularly, using nu- 
merical simulations we determine the quantum phase diagram 
of a hardcore (HC) bosonic t-Jz planar (2D) model, illustrat- 
ing the main conclusions of this paper. Our control parameter 
is the lattice anisotropy 7, < 7 < 1, with 7=1 representing 
the isotropic case with no explicitly broken lattice rotational 
symmetry 0]. In the small 7 regime a confinement interac- 
tion emerges with a non-vanishing topological hidden order 
and static incommensurate magnetic orders. This is the static 
stripe phase which, as the anisotropy is increased, persists up 
to a critical value 7c where a confinement-deconfinement tran- 
sition occurs. At that point the topological stripe order melts 
and gives way to a phase-segregated state for > J^^ (for 



7 = 1 see 01) or to a superfluid phase with static AF correla- 
tions and dynamical stripe fluctuations for > > J^^ . 

There is a number of significant reasons to be interested in 
this model. First of all, its fermionic counterpart is believed 
to be a relevant model for electron motion in the copper oxide 
planes of HTSCs, responsible for unique properties of these 
materials. Second, we analytically prove that the low-energy 
spectra of fermionic and bosonic versions of the model are 
identical in the limit of weak coupling 7 between the chains. 
This is a natural consequence of an exact result in one dimen- 
sion 19]. Finally, due to the absence of the infamous /erm;o« 
sign problem, we are able to simulate the properties of the 
model for large lattice sizes and very low temperatures. 

We consider a 2D anisotropic t-Jz model on a square lattice 
for the spin-1/2 HC bosons defined in Ref. 0] 



E^'feWa+H.c.) , (1) 
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where vectors i = x ex + y ey run over all the Ns sites of a 2D 
square lattice, e^, are the unit vectors of this lattice, and bars 
over the operators imply that the constraint of no more than 
one particle per site (the HC constraint) has been already in- 
corporated into the operator algebra. The spin and the number 
operators on site i are defined by the following expressions 



nil) and n; = + rii^ 



(2) 



with riia- ~ ^L^icr- constrained HC boson operators 

can be expressed in terms of standard HC boson operators, 
foL, which are HC in each flavor and obey the com- 



mutation relations: 



[bl,bl,] ^ 0,[b,^,bl,] = 
The connection between the two bosonic 
algebras is given by = ^[^.(l — Uig), where a denotes 
the spin orientation opposite to a. From this expression we 
can derive the commutators for the algebra which defines our 
spin-1/2 HC bosons 
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Using the general transformations introduced in Refs. 
we can map these particles into the HC-constrained fermion 
operators, c|„ = cj^(l - nnf) and = (1 - "-iajq^ 



cxp[ 



(4) 



Here, is the angle between the spatial vector j — i and 

a fixed direction on the lattice. Given this transformation, Hb 
may be rewritten in terms of the fermionic operators 



Hp 



H.c. 
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withyl^(i) = ^j,[Lj(k, i) — tj(k, i + e^)] nj.. Expressions Q 
and (|5ji for Hb only differ in the kinetic-energy term, with the 
fermionic language including a non-local gauge field A^{i) 
associated with the change in particle exchange statistics. In 
general, this gauge field cannot be eliminated in dimensions 
larger than one which means that different particle statistics 
can give rise to different physics. However, we will see be- 
low that the gauge field is irrelevant and can be eliminated 
for the low energy spectrum of Hb in the strongly anisotropic 
region It**! <C \t'^\- Notice that can be finite and there- 
fore the problem is still 2D. In other words, when the ratio 
7 = l^^l/l^^ l is small the properties of the bosonic model 
governed by the low-energy spectrum of Hb will be identical 
to the corresponding properties of the fermionic t-J^ model 
[Hp = HeiA^ii) = 0)]. 

Let us start by considering the limit = 0. In this limit, 
the system is an array of one-dimensional (ID) t-J^ chains 
which are only magnetically coupled through J|. From the 
quasi-exact solution of the ID i- model it is known that 
the lowest energy subspace Aio consists of states in which 
the spins are antiferromagnetically ordered and each hole is 
an antiphase domain for the Neel OP (see Fig. la). Since the 
excited subspaces containing an extra spin excitation are sep- 
arated by a finite energy gap, Mq becomes the relevant sub- 
space to describe the low-energy physics of weakly coupled 
chains. The combination of the interchain magnetic coupling 
and the fact that each charge is an antiphase domain gives rise 
to a confining interaction between holes on adjacent chains. 
This is because the misalignment of holes on the neighboring 
chains either breaks the AF bonds or introduces ferromagnetic 
links, both energetically unfavorable (Fig. QJ)). The slope of 
the linear potential is proportional to = 7 . This confin- 
ing interaction leads to a formation of a hole stripe, which is 
simultaneously a 2D antiphase domain boundary for the Neel 
OP fioll . Note that the characteristic length of stripe fluctua- 
tions is ^ ~ 7~^/'^. For finite hole densities, ph, the stability 
condition for the existence of a stripe phase is 7 > t^/ Jf^. 

What happens when a small hopping <C , is in- 
cluded? As shown in Fig.^, the hopping of the hole to the 
adjacent chain has an energy cost of /2. Since <gi J^, 
the magnetic structure again acts as a potential barrier confin- 
ing the hole to move in the a;-direction. The main effect of 
is a second-order diagonal correction, 2(i^)^/ (see Fig. Ic), 
to the energy of the hole. In other words, the low-energy ef- 
fective model for <^ , is essentially the same as the 



one for = Q. Consequently, the low-energy spectra of the 
bosonic and the fermionic t-J^ Hamiltonians are the same for 
ty <C t^, J2 , i.e., the gauge field of Eq. (|5jl can be eliminated 
at low energies because the hole is effectively moving in the 
x-direction (the transmutator of statistics is a symmetry 
of Hb restricted to its low-energy sector). This concludes the 
proof of the equivalence of the low-energy spectra of bosonic 
and fermionic 2D t-J^ models in the strongly anisotropic case. 



f (2|a:| - 1 + (S(i)) 




(!>)■ 



FIG. 1: Illustration of the intrinsic stripe formation mechanism in 
the 2D anisotropic t-Jz model (a), and energy costs associated with 
different processes: intra-chain hopping (b) and inter-chain hopping 
(c). Dynamical confinement emerges from the competition between 
magnetism and delocalization. 

Since a model defined by the Hamiltonian Q is bosonic, its 
quantum Monte Carlo (QMC) simulation is not affected by the 
sigyi problem. We have used the worldline loop algorithm in 
the continuous imaginary-time formulation (for a review see 
iHHi ). similar to the one used for the t-J model in Ref. fl2il . 
This efficient method allows one to perform calculations at 
very low temperatures and large lattice sizes. Unless noted, all 
simulations were performed on square lattices of linear size L 
ranging from 10 to 40, with ph = 0.2, and error bars smaller 
than the size of the symbols. The anisotropy was parame- 
terized by setting \t^\ = 1 (all other energy parameters are 
measured in units of \t^'\) and introducing a single parameter 
7 — \ty\ = jy/J^, which may vary between zero (discon- 
nected chains) and unity (isotropic 2D model 0]). Tempera- 
ture was fixed by putting /3 = |P|/(A:bT) = 50, which was 
low enough to sample essentially the properties of the ground 
state. The method is exact up to a statistical error and we are 
using periodic boundary conditions to avoid the artificial os- 
cillations introduced by open boundaries Iil3»1 . 

What observable quantities characterize a stripe phase? As 
argued above one needs to measure both the magnetic struc- 
ture factor (MSF) and a topological hidden order to uniquely 
determine it. The MSF {N^ = L^) 
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is expected to display incommensurate peaks at wavevectors 
k = Q ± (ttJ, 0) in the presence of stripes oriented along 
the y-direction (across the chains), reflecting the 2D character 
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of the new antiphase boundaries (see Fig. ^). Here Q = 
(tt, tt) is the AF wavevector The results for the AF OP are 
presented in Fig.|2lfor = 1. For small values of 7, 7 = 
0.2 and 7 = 0.4, there are pronounced peaks at wavevectors 
k = (tt ± 7r/5, tt) indicating a superstructure with a period 
of ten lattice spacings, consistent with incommensurate spin 
ordering and a stripe-like phase. At 7 = 0.6 the height of the 
incommensurate peaks drops while, simultaneously, the peak 
at k = Q starts to grow, and becomes dominant for 7 = 0.8, 
signaling a commensurate AF spin ordering. 
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FIG. 2: MSF S{k) for a 30 x 30 system at hole doping p,, = 0.2, 
— 1, and different values of anisotropy 7. Error bars have been 
omitted for clarity, standard relative error for any data point never 
exceeds 10%. 

The presence of incommensurate peaks in S'(k) does not 
unequivocally prove the existence of a stripe phase since a 
similar behavior can be obtained, for example, for a spin spiral 
phase. It is therefore crucial for our analysis to define another 
quantity which unambiguously signals the presence of stripes. 
Since stripes are topological defects (antiphase boundaries), 
this quantity is called topological OP (TOP) fl4!] and is calcu- 
lated for every ID chain with fixed y coordinate. The corre- 
sponding two-point correlation function is 



Giy) 
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where S^^ is the spin projection operator at site i = a; e^; + 
y e.y, and Np = L{1 — ph) is the number of particles in 
the chain. In a system with periodic boundary conditions 
Eq. is, of course, independent of y, allowing us to av- 
erage the results of TOP measurements for different chains. 
The parameter T{x,x',y) is defined by T{x,x',y) = x' + 

12p=o ^ ^{x+p,y))- Without the second term this param- 
eter would just turn Eq. Q into the square of the ID Neel 
OP. The second term, however, introduces an additional fac- 
tor of (—1) for every hole encountered between the positions 
X and x' in the chain with fixed coordinate y, indicating that 
the hole is an antiphase boundary for the Neel OP. It can be 
easily established that the TOP defined by Q will reach its 



maximum value of unity when evaluated in the ground state 
of the ID t-Jz model. Thus, G{y) quantitatively measures to 
what extent the separate chains in the 2D system have retained 
their characteristic ID ground state topological features. A 
non-vanishing TOP and incommensurate spin ordering in the 
MSF constitute strong evidence that the system is in the stripe 
phase. 

In order to rule out the possibility that the observed behav- 
ior is a finite-size effect, we have extrapolated the values of 
the MSF at the position of the incommensurate peak (Fig.|3 
and TOP (Fig.0} for small values of 7 to the infinite system 
size. Both quantities clearly tend to finite values in the ther- 
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FIG. 3: Finite-size scaling of the MSF 5(k) at the position of the 
incommensurate peak k — (47r/5, tt) for = 1 and different val- 
ues of the anisotropy 7: 0.1 (circles), 0.2 (squares), 0.3 (diamonds), 
0.4 (triangles). Lines are linear fits to the QMC data. Data points in 
the inset are obtained by extrapolation to the L ^ 00 limit, with the 
dashed line used as a guide to the eye. 

modynamic, L ^ 00, limit for the range of 7 studied. It is 
noteworthy, that for the same 7 range the value of the MSF 
at the AF wavevector Q extrapolates reliably to zero (within 
error bars). 

The non-uniform behavior of the incommensurate peak 
height (inset of Fig.|3} is a consequence of two competing pro- 
cesses. A stronger coupling between chains leads to a stronger 
confining potential, facilitating the formation of stripes. On 
the other hand, the same increase tends to disturb the ID 
ground state of the chains due to inter-chain hopping and in- 
teractions, driving the system away from stripe ordering. This 
competition leads to an optimum value 70 ^ 0.3, for which 
the peak height reaches its maximum. 

At 7 = 7c (see Fig.|5jl, both the TOP and the incommen- 
surate peak in S'(k) extrapolate to zero in the thermodynamic 
limit. For 7 > 7c, we observe nonzero values of the super- 
fluid density ps with dominant AF correlations if > J^^ 
(for Jz < J^^ the system is in the Nagaoka ferromagnetic 
superfluid state), ps was calculated in a 10 x 10 cluster by 
computing the mean square deviation of the spatial winding 
number 1 15]. Coexistence of AF correlations and superfluid- 
ity indicates that the bosons could be moving in pairs since 
the propagation of individual bosons destroys the AF order- 
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FIG. 4: Finite-size scaling of the TOP G for = 1 and differ- 
ent values of the anisotropy 7: 0.1 (circles), 0.2 (squares), 0.3 (dia- 
monds) and 0.4 (triangles). Lines are second-order polynomial fits to 
the QMC data. Data points in the inset are obtained by extrapolation 
to the L — > 00 limit, with the dashed line used as a guide to the eye. 

ing. For > ~ 0.5, there is a quantum phase transition 
to a phase-segregated state where a hole-rich superfluid phase 
coexists with an AF region with no holes. This is the origin of 
the sharp AF peak of Fig.|2l(7 = 0.8). Most likely, the glue 
that keeps particles together in the hole-rich region is still the 
confining interaction provided by the magnetic background. 
In summary, we have determined the quantum phase dia- 



gram of the planar anisotropic bosonic t-J^ model with peri- 
odic boundary conditions. For small values of the anisotropy 
parameter 7 and hole doping ph, there is a quantum stripe 
phase characterized by topological hidden order, and incom- 
mensurate peaks at k = Q ± {tt6, 0) in the MSF The stripes 
are formed in the absence of any physical long-range interac- 
tions. The glue stabilizing each stripe is a confining poten- 
tial that emerges dynamically out of the competition between 
the kinetic energy and the AF exchange. For anisotropics 
7 > 7c a quantum phase transition to a superfluid (decon- 
fining) phase, coexisting with AF coiTelations, takes place if 
J^^ > > J^^. For larger values of the system phase 
segregates. It is remarkable that the same attractive interac- 
tion gives rise to both the stripe phase and phase segrega- 
tion. The existence of a transition between a superfluid and 
a phase-separated state in proximity of the stripe instability is 
very suggestive when compared with the phenomenology of 
the HTSC. The phase-separated region must survive under a 
transmutation of the statistics because the fermionic kinetic 
energy is always higher than the bosonic one. The presence 
of an attractive interaction that segregates the fermions can 
easily induce a superconducting state if is lower but close 
to the critical value that leads to the phase-segregated state. 
Note that for the bosonic case, this critical value {J^^ ~ 0.5) 
is slightly larger than the values of that are considered re- 
alistic for the cuprate HTSC. 

J.S. is grateful to the Swedish Foundation for Strategic Re- 
search, Swedish Research Council, Goran Gustafsson Foun- 
dation, SNAC, and SHARCNET. 



[1] J. Zaanen and O. Gunnarsson, Phys. Rev. B 40, 7391 (1989); K. 

Machida, Physica C 158, 192 (1989). 
[2] J.M. Tranquada etal.. Nature 375, 561 (1995). 
[3] C. D. Batista, G. Ortiz and A. V. Balatsky, Phys. Rev. B 64, 

172508 (2001). 

[4] B. Normand and A. R Kampf, Phys. Rev. B 64, 024521 (2001) 

studied fermionic Hamiltonians for 7 ~ 1. 
[5] M. Boninsegni, Phys. Rev. Lett. 87, 087201 (2001). 
[6] C. D. Batista and G. Ortiz, Phys. Rev. Lett. 85, 4755 (2000). 
[7] C. D. Batista, G. Ortiz and J. E. Gubernatis, Phys. Rev. B 65, 

180402 (2002). 



[8] C. D. Batista and G. Ortiz, Phys. Rev. Lett. 86, 1082 (2001). 

[9] C. D. Batista and G. Ortiz, Adv. in Phys. 53, 1 (2004). 
[10] J. Zaanen etal., Philos. Mag. 81, 1485 (2001). 
[11] H. G. Evertz, Adv. Phys. 52, 1 (2003). 
[12] B. Ammon et ai, Phys. Rev. B 58, 4304 (1998). 
[13] S. R. White and D. J. Scalapino, Phys. Rev. Lett. 80, 1272 
(1998). 

[14] M. den Nijs and K. Rommelse , Phys. Rev. B 40, 4709 (1989). 
[15] E. L. Pollock and D. M. Ceperley, Phys. Rev. B 36, 8343 
(1987); S. Zhang etal., Phys. Rev. Lett. 74, 1500 (1995). 



4 



LA-UR-03-XXXX February 2, 2008 submitted to Physical Review Letters 




7 

FIG. 5: Quantum phase diagram of the anisotropic bosonic t-J^ 
model for ph = 0.2 (p| 7 1)' displaying pictorial represen- 
tations of the phases. 
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